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Abstract. Using an Environmentally Friendly Renormalization Group we derive an 
ab initio universal scaling form for the equation of state for the O(N) model, y = f{x), 
that exhibits all required analyticity properties in the limits x — > 0, x — > oo and 
x — > — 1. Unlike current methodologies based on a phenomenological scaling ansatz 
the scaling function is derived solely from the underlying Landau-Ginzburg- Wilson 
Hamiltonian and depends only on the three Wilson functions 7a, j v and 7^2 which 
exhibit a non-trivial crossover between the Wilson-Fisher fixed point and the strong 
coupling fixed point associated with the Goldstone modes on the coexistence curve. 
We give explicit results for A = 2, 3 and 4 to one- loop order and compare with known 
results. 



PACS numbers: 64.60.Ak, 
1. Introduction 

The equation of state for the O(N) model has been an object of intense scrutiny over the 
last 30-40 years (see, for instance, [1] and [2] for recent reviews). It exhibits crossover 
behaviour between three distinct asymptotic regimes - the critical region approached 
along the critical isotherm, the critical region approached along the critical isochore, 
and, finally, the coexistence curve. It is the difficulty of encapsulating these distinct 
scaling behaviours within one overall scaling function that has prevented its ab initio 
derivation from an underlying microscopic model. 

The first attempts at such an ab initio calculation, using the renormalization group 
(RG) and an e-expansion [3], foundered on the fact that they did not exhibit Griffiths 
analyticity in the large x = t/tp 1 ^ limit. This was due to the fact that the expansion was 
around a particular fixed point - the Wilson-Fisher fixed point. However, to obtain a 
universal equation of state, valid in the entire phase diagram, using the RG and without 
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phenomenological input, the 93-dependent crossover between this fixed point and the 
strong coupling fixed point associated with the massless Goldstone excitations on the 
coexistence curve must be accessed and controlled. In the latter regime the non-linear a 
model gives a good description [5] . However, this model only accounts for the Goldstone 
bosons and does not offer a full description of the phase diagram. 

As pure first principles calculations using a Landau- Ginzburg- Wilson Hamiltonian 
and the RG have not been capable of obtaining an equation of state valid in all 
asymptotic regimes, the current "state of the art" [U E] is to base calculations on 
a parametrised phenomenological scaling form [6] that has the correct analyticity 
properties. Instead of the magnetisation, <p, reduced temperature, t, and magnetic 
field, H, new variables, 9 and R, are introduced, the relation between them being 



The particular functional dependence on R ensures that Griffith's analyticity is 
preserved. However, the two functions m{9) and h(9) are arbitrary. Hence, any 
particular choice constitutes a pure ansatz. Thus, Schofield's scaling form gives a large 
class of models, all of which, by construction, are consistent with Griffith's analyticity. 
How a particular microscopic model is related to this large class depends on how these 
undetermined functions are represented, the canonical approach being to represent them 
as polynomials. Thus, there is no unique map between a microscopic model and a 
member of the class of Schofield parametrisations. 

In this new parametrisation the scaling function of the equation of state is given by 



Most current field theoretic formulations for determining the equation of state (see 
[2] for a comprehensive review) rely on such formulations. The drawback is that the 
underlying microscopic theory is not used to determine the functional form of m(9) and 
h(9). Rather, an ansatz is made as to the general functional form, which depends on 
certain unknown parameters, and then the underlying microscopic theory is used to 
fix them. The most common ansatz is that the functions are polynomials in 9. The 
coefficients of the powers of 9 are then determined by calculating the values of certain 
observables independently from the underlying microscopic theory. 

In contrast to the above, in this papejj], using only the Landau- Ginzburg- Wilson 
Hamiltonian for the O(N) model and by implementing an Environmentally Friendly 
RG [8j [9] that tracks the crossover between the fixed points that control the different 
asymptotic regimes, we derive an universal equation of state that obeys all required 
analyticity properties and where no phenomenological input is required, only the three 
Wilson functions 7^, 7^2 and 7a- 

f Which is a follow up to [10] which treated the case N = 1 







where 




(1) 
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2. The Equation of State 

In the critical region, the equation of state [H] relates the external magnetic field H, 
the reduced temperature t, and the magnetization (p, 

y = f(x) (2) 

where the universal scaling function f(x), is normalised such that /(0) = 1 on the 
critical isotherm, and /(— 1) = on the coexistence curve. The scaling variables y and 
x are given by y — B^H/tp 5 and x = B 1 '^t/ip> 1 ^ . 

The function f(x) has expansions around the limits x = and x —>■ oo given by 



CO 



fix) = 1 + E fn* n f(^) = ^E fnX- 2n " (3) 

n=l n=0 

In the limit x — > oo a natural variable is 2 = biip/t^, where &i an amplitude ratio in 
terms of which the equation of state takes the form H oc t l3S F(z), where the universal 
scaling function F(z) for small z has an expansion of the form 

F M = ' + i* ,+ Scs£i)i 1 ~ (4) 

where r2 = r4 = 1 by choice of normalisation. As (jlj) is an expansion in (p, the 
constants r 2n are related to the 2n-point correlation functions at (p — and hence 
are very natural observables to calculate in lattice simulations. In the limit z — > oo, 
F(z) has an expansion of the form 

CO 

F(z) = z s E Fr^ fc//3 (5) 

fc=0 

By relating f(x) and i^(^) expansion coefficients of the two functions can be related to 
find 

foo _ 2n+l-£ r 2ra+2 / fi s 

Jn ~ F °°(2n+1)! 1 j 

poo 

f _ f_n_ -n/P / ? \ 
/n — P co Z I' J 

^0 

where is the universal zero of the equation of state in terms of the variable z. Thus, 
we see it is sufficient to know the expansion coefficients of f(x) in the limits x — ► and 
oo in order to calculate the asymptotic properties of F(z) and the interesting functions 

T2n- 

Unlike the limits x — > and oo, near the coexistence curve, x — > —1, there are 
no rigorous mathematical arguments as to the analyticity properties of f(x), although 
there do exist conjectures. In [12], based on an e-expansion analysis, it was conjectured 
that (1 + x) has a double expansion in powers of y and y( d_2 )/ 2 of the form 

1 + x = Cl y + C2 y 1 ~ e/2 + d lV 2 + d 2 y 2 ~ e/2 + ... (8) 

In three dimensions it predicts an expansion of (1 + x) in powers of y 1 ^ 2 . Studies of the 
non-linear a model lead one to expect a leading behavior of the form 

f(x) ~ c,(l + xf/^ (9) 
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though, as mentioned, the nature of the corrections to this behavior is not well 
understood, although (jSJ) is one conjecture. In the 1/N expansion there is some evidence 
[2] for logarithmic corrections of the form ln(l + x) in three dimensions. 

3. A Renormalization Group Representation of the Equation of State 

In this section we briefly outline the derivation of the equation of state for a theory 
described by the standard Landau- Ginzburg- Wilson Hamiltonian with O(N) symmetry 

H[ip] = J d d x Qv^V^+ ^r(*W+ ^>V) 2 ) (10) 

with r = r c + t B , where r c is the value of r at the critical temperature T c and 
t B = A 2 ( - T ~ Tc ^ , A being the microscopic scale. Due to the Ward identities of the 

J- c 

model, it is sufficient to know only the transversal correlation functions y\ N ' M \ as all 
the other vertex functions can be reconstructed from these. For instance, the equation 
of state itself is given by H = T t ip. 

Due to the existence of large ffuctations in the critical regime a renormalization of 
the microscopic bare parameters of the form 

t(mt, «) = Z~£{K)t B {mt) (11) 
A(k) = Z x {k)X b (12) 
<p( K ) = Z- l l\K)ip B (13) 

must be imposed, where k is an arbitrary renormalization scale and m t is the inverse 
transverse correlation length. The renormalized parameters satisfy the differential 
equations 

K dtK /.^^ where = ^.Z^\ c (14) 

an uk 

K ~dK = ^( K )K K ) where jx(k) = k-^- In Z x \ c (15) 

d<f(K) l/\/\ d . 

K — = - ■^J'pWVW where "f v (K) — k— in (16) 

where the right hand side are the Wilson functions associated with this coordinate 
transformation and the derivative is taken along an appropriately chosen curve in the 
phase diagram, which we here denote by c. 

Integration of the RG equation for any multiplicatively renormalizable y[ n ' m ^ yields 

r\ N ' M \t, a, cp) = jn^rr^), A(«), (17) 

The renormalization constants Z v , Z^ and Z\ are fixed by imposing the explicitly 
magnetization dependent normalization conditions on the transverse correlation 
functions 

dp^t{p, «), A(/c), k) 

^(O,*^^)^^), k) =1 (19) 

r< 4) (0, *(«,«), A(«), =A. (20) 



p 2 =0 ^ 
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Figure 1. Integrating along curves of constant (p (or t) in the phase diagram 
starting from the coexistence curve. 



while the condition 



K 2 



rf (0,*(«,k),A(k), p(k),k), (21) 



serves as a gauge fixing condition that relates the sliding renormalization scale k to the 
physical temperature t and magnetization ip. Physically, k is a fiducial value of the 
non-linear scaling field m t . 

Besides m t , the other non-linear scaling field we use to parametrise our results is 

2 i rSV 

™ v = -— w (22) 

6 o P ^t V=o 

which is a RG invariant. It represents the anisotropy in the masses of the longitudinal 

j p 2l t \p 2 



and transverse modes and is related to the stiffness constant p s = <p 2 d p 2 r| ^ |p2_Q via 



m 



2 



\Xp s - With this renormalization prescription one may determine the equation 
of state in terms of the non-linear scaling fields m t and the transverse and 

longitudinal propagators that appear in all perturbative diagrams can be parametrised 
in terms of them. 

We now give a short review of the general methodology used to derive the equation 
of state as developed in reference 0, HD]- The equation can be found by integrating 
dTt(t,(p) along a curve of constant ip (see figure ([!])), where in dt = dV[ 2 ' /I"^ 2 , the 
right hand side is written in the coordinate system (mt,m v ). We note the following 
points [TO] : 

1) We take m f = k as the flow variable of the RG and hold ip constant. The 
non-linear scaling variable, m v , is then k dependent through its definition ( |22l) . 

2) We integrate dt along a curve of constant ip and fix the boundary condition on 
the coexistence curve, where m t = 0, whence we may write (T — T c (ip)) = t + A, where 
A = (T c — T c (ip)) is the temperature shift that measures the distance between the critical 
point and the point on the coexistence curve, T c (ip), see figure ([TJ. One finds 

A x (l + x) = F(z) (23) 
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where z = mt/m^ and the scaling variable x = B l /^t/ip 1 /^ , B being a amplitude. The 
universal scaling function T{z) is given by 

2- 7A +7 ¥ > y idx , . 

w 7o 2 - 7A + 7v x y ' 

where we have defined Aji = — 7^ F ), WF denoting the Wilson-Fisher fixed point. 

3) The equation of state is obtained using H = T^(p and the RG equation for . 
One obtains H/ip s = A^ 1 Q(z), with the universal function 

7 C z A 7A~ A ~>V> dy _ P* 2&~ft£ dy 

Q{z) =qI 3 J^> 2 ~~t\+~tV VQ Joo y (25) 

Defining the scaling variable y = A 3 1 (H/ip 5 ) we then obtain the universal equation of 
state in the form ([2]) with 

f{x) = ^g{F-\A l {l + x)). (26) 

4) The amplitudes A\ and A3 are both universal and can be determined as 
A\ = J-(z c ) and A 3 = G(z c ), where z c is the value of z that corresponds to the critical 
isotherm. 

In order to determine the expansion coefficients /° and f£°, as introduced in section 
1, one requires the Taylor expansion of f(x) around x = and x = 00. In terms of 
our parametric representation, d n f(x)/dx n can be expressed using d/dx = (dz/dx)d/dz, 
where dz/dx = Ai/(dJ r (z)/dz) hence, 

dPf{x) A 1 (fd^{z)\- 1 dY 

-^-T 3 {[-ir) Tz) G{z) (27) 

which needs to be evaluated at the points of interest z = z c (x = 0), z = 00 (x = 00) 
and z = (x = — 1). 

4. Results 

4-1. The one loop Wilson functions 

To one loop order the bare vertex functions are given by 

if = p > +r + ^ + ^ < 3 > +±(N-1)Q, (28) 

r ! 2 " = 1 - 1 (O + ) < 29 > 

if = A - §A> ( + ) (30) 

where the dashed line in the loops denotes the transversal propagator and the continuous 
line the longitudinal one, i.e., /~\ = f -M-4? } x and /~\ = f Jr4x , rm 
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Using (TT8|) . f)19p and (1201) . the renormalization constants can be determined and hence 
the Wilson functions 

d ( N 



1* =0 (32) 
A d ( N - 1 \ 

The running dimensionless coupling A satisfies 

d\(z) , l2/ , A 1 ,d-a (iV — 1)\ , . 

Taking the initial condition A(zo) — A, in the limit zq — > oo, A — > oo one arrives at the 
universal separatrix solution^] 



A W = [^[(1 + 1)^ + ^1 ] (35) 



On the separatrix 



1^ i (N-l) 



(1 + ?)^ + 



a =(| ' (/)| (36) 

7^ =0 (38) 

With the Wilson functions in hand we can calculate the scaling functions. In three 
dimensions the scaling function Q(z) is analytic, while the function T{z) can be written 
as an integral. Explicitly, 



3 / /l + l/ z 2_ \ TO 

GW=**Il r x _ ri (39) 



1=1 



12(iV-l) 



T{z) = f 
Jo 



'l + v /l + l/ x 2\ n^tojw+ii 



o 2 -7 A 



(/ \ n i 

M + 1/ X 2_ r \ 
± I x—dx (40) 

where r« satisfies the cubic equation (N — l)rf + 18rf — 9 = and rii = 

9(l-N)+(161+2N-N 2 )ri+9(N-l)r^ 4 

I2r t + (N-I)rf (10-N)(N+8) " 

In the limit z — > oo, which corresponds to approaching the critical point along the 
critical isochore, the Wilson-Fisher fixed point is approached and 7$ — > 7,^ with, at 
one loop, 7a = (4 — d) and 7^2 = (4 — d)(N + 2)/(N + 8). On the contrary, in the 
limit z — > 0, which corresponds to approaching the critical point along the coexistence 

§ This solution may also be reached by choosing the initial coupling to be on the separatrix solution 
at z = z . 
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curve x — > — 1, the strong-coupling fixed point is approached and 7$ — > 7? . There, 
the Goldstone bosons dominate and 7a = 7^,2 = (4 — d). Finally, the critical isotherm, 
x = 0, is reached in the limit z — > z c . 

With (jlQjl and (139|) we can plot the scaling function, /(x), for the full universal 
equation of state, as seen in Figures El E] and H] for the cases N = 2, 3 and 4 respectively. 




-1 -0.5 0.5 1 1.5 



Figure 2. The scaling function f(x) for the three-dimensional XY model. The 
dashed line was taken from reference [15j . 




-1 -0.5 0.5 1 1.5 2 



Figure 3. The scaling function f(x) for the three-dimensional Heisenberg 
model. The dashed line was taken from reference [13j . 



4-2. The limit z — ► 00 

In the limit z — > 00, the Wilson functions can be expanded as power series in z~ 2 

00 

7,(z)=7r + E^K 2n (41) 

71=1 




Figure 4. The scaling function f(x) for the three-dimensional 0(4) model. The 
dashed line was taken from reference 



Hence, the universal scaling functions JF and Q can also be written as power series in 
z~ 2 . This is true in a diagrammatic expansion to all orders, not just at one loop. In the 
limit z — y 0, ji — > but the nature of the corrections is not obvious. At the one-loop 
level, from (!36|) . one can see that the leading corrections to the strong coupling fixed 
point values will be z^~ d ^ 2 . 

To determine the constant A\ we take the z — > oo limit of (|40p . identify the divergent 
part with A\X and the constant remainder with A\ ijj] For instance, for N = 1, 2, 3, 4 we 
can write the divergent component of the function T: 

J ~= 4 - 4 Hs^r)? (42) 

2 " + { 4 +10 t ± r.in -l)(r i + 12) i ~* 



r - 4 + 2 + 18-158r,-18rM J_ 

3 + V 77 + 77 ^ r< (r 4 - l)(2r 4 + 12) J *» ^ 

C = 4 + f-Z + ± f 27-153r,-27rn J_ 

4 { 2 + 77^r J (r J -l)(3r J + 12)J *> ^ 

The universal amplitudes v4 3 can then be obtained using the corresponding equation 
in 4) in section [3j With these two constants, along with the two scaling functions T 
and Q the scaling function of the universal equation of state may be obtained, and, 
correspondingly, any expansion coefficient, /°, or r n . 

In Table [TJ we see some results for a variety of expansion coefficients associated 
with the asymptotic regimes x — > and x — > oo. We also compare with various 
other known results. The agreement is as good as one would expect from a one-loop 

|| If 2(4 — d)/3(d — 2) is a positive integer, n, then this remainder is zero and A\ cannot be determined 
by looking at the asymptotic limit z — > oo. This is a pure artefact of the one-loop approximation, 
where r\ = 0, and has no physical meaning. 
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approximation. However, of note here is not the precision of the estimates but that 
they have been obtained from an ab initio calculation where no phenomenological input 
was necessary. Thus, any expansion coefficient is obtained by an expansion in the 
appropriate asymptotic limit of the universal functions T{z) and Q{z) which in their 
turn depend only on the Wilson functions. 

Table 1. Values of some one- loop amplitudes for the 3-dimensional O(N) 
model. We have written some known results as reported in [2] (see table (21)) 
for N=2; [13j for N=3 and [14] for N=4. 



iV 


o 
i 


Q 

o 


4 


Zc 


0.694 


0.677 


0.656 


At 


0.959 


0.934 


0.897 


A 3 


1.022 


0.840 


0.698 


R\ 


6.35 ( 7.6(2) ) 


6.06 (7.8(3) ) 


8.86 ( 7.6(4) ) 


ipoo 


0.035 ( 0.0304(3) ) 


0.0348 (0.0266(5)) 


0.015 ( 0.0240(5) ) 


R x 


1.42 (1.41(6)) 


1.280 (1.31(7)) 


1.18 ( 1.12(11) ) 


Jl 


0.746 


0.789 


1.246 


foo 
J2 


0.641 


0.695 


1.171 


foo 
J3 


0.103 


0.097 


0.151 


foo 
JA 


-0.005 


-0.006 


-0.014 


C f 


15.6 (15(10)) 


4.92 ( 5(3)) 


2.66 ( 2.8(1.4) ) 


A° 


1.16 


1.26 (1.34(5)) 


1.35 ( 1.5(8)) 


fS 


0.086 


0.143 (0.020(2)) 


0.2031 ( 0.33(5)) 


JS 


-0.022 


-0.040(-0.10(l)) 


-0.056 (-0.08(2) ) 


re 


2.704 ( 1.951(14) ) 


2.90 ( 2.1(6)) 


2.11 ( 1.79(2) ) 


r» 


2.873 ( 1.36(9) ) 


2.80 ( 0.6(2) ) 


1.29 ( 0.2(4)) 


no 


-1.58 ( -7(5) ) 


- 2.05 ( -6(3) ) 


-0.97 ( -5(6)) 



4-3. The limit z —>■ 

In the limit z — > 0, corresponding to the coexistence curve, the parametric functions 
T{z) and Q{z) can be written as power series in the scaling variable z as T{z) = 
A X {1 + x) = %z 2 + \z z + f z 4 + . . . and Q{z) = A 3 y = alz 4 + fz 5 +gz 6 + .... We can 
invert the equation for T{z) to obtain the non-linear scaling field z in terms of (1 + x) 
and then, once we substitute into the function Q(z), we obtain the equation of state 
near the coexistence curve in the form 

y = c f (l + xf + ci(l + x) 5/2 + c 2 (l + x) + ... (46) 

or alternatively 

(1 + x) = c x y + ~c 2 y l/2 + c 3 y 3 ^ + ... (47) 
In Table [2] we give the first expansion coefficients and compare them to the results of 
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Table 2. Values of some one-loop amplitudes for the 3-dimensional O(N) model 
in the limit z — > 0. The values in parenthesis for c\ y C2 were taken from [12j. 
where an e expansion was used. 



N 


2 


3 


4 


Cl 


-111.6 


-17.8 


-6.4 


C2 


399.8 


32.0 


7.6 


Cl 


0.82(0.9) 


0.66(0.82) 


0.54(0.75) 


C2 


0.25(0.1) 


0.45(0.18) 


0.61(0.25) 


C3 


0.46 


0.55 


0.57 



5. Conclusions 

Although there exist established methods for calculating the universal equation of state, 
f(x), while preserving Griffiths analyticity, they are based on a phenomenological scaling 
ansatz that has no underlying microscopic basis. On the other hand it has not been 
possible to preserve all asymptotic properties of f(x) starting from a Landau- Ginzburg- 
Wilson Hamiltonian using the RG, due to the fact that the latter involved an expansion 
around the Wilson-Fisher fixed point. In order to calculate f(x) using RG methods it 
is necessary to be able to capture dependent crossover between the two different fixed 
points. 

In this paper we have developed a RG that captures the crossover between the 
Wilson-Fisher and strong-coupling fixed points and therefore captures all elements of 
Griffiths analyticity, accessing all three different scaling regimes. Thus, we have achieved 
an ab initio calculation of the universal equation of state from an underlying microscopic 
model with no phenomenological input. The calculation depends only on the three 
Wilson functions 7^, 7^,2 and 7a- We used the method to calculate f(x) to one loop 
for arbitrary N and made a comparison to known results from other methods. The 
novelty of the present approach is not in the precision of any result for an expansion 
coefficient, as here we have only worked to one loop, but rather in exhibiting a method 
that order by order in perturbation theory preserves all necessary analyticity properties 
of f{x). However, the method can be extended to higher loop order. The chief difficulty 
in doing so is that the higher order Feynman diagrams must be evaluated numerically 
(remembering that they are crossover functions not constants) as there are no closed 
form expressions for them. As the amplitude Ai involves cancelling two divergent 
expressions this also is trickier when done numerically. We will return to the question 
of two-loop calculations in a future publication. 
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